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АНАЛИЗ МАТЕМАТИЧЕСКОЙ МОДЕЛИ КОЛЬЦЕВОЙ РОТОРНОЙ ПЕЧИ 


Описана дискретная модель температурных полей кольцевой роторной печи. Составлена компьютерная про- 
грамма для рас чёта тепловых процессов и построена передаточная матрица рассматриваемого объекта. 
С помощью спект; ров Г ершгорина исследованы характеристики передаточной матрицы тепловых полей рас- 
сматриваемой печи как объекта управления. 
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Введение. В современных условиях наблюдается острая нехватка строительных материалов. 
Разработка нового современного высокопроизводительного оборудования для производства 
строительных материалов, безусловно, является актуальной задачей. Для производства извести 
применяют печи различных типов: шахтные, вращающиеся и др. Кроме этого используют 
установки для обжига извести во взвешенном состоянии, в кипящем слое, на специальных решёт- 
ках ит.д. 

Кольцевая роторная печь (рис.1), как показывают расчёты, является более энергосбере- 
гающей по сравнению с традиционными печами, применяемыми для производства 1 т извести. 
Основная задача технологического процесса заключается в термической обработке исходного 
материала (щебня), перемещаемого вдоль печи специальным устройством. 














ых 








Рис.1. Схема объекта управления: 1 - вход и выход камеры; 2 — форсунки; 3 -— датчики, 
измеряющие температуру в камере; 4 -— кирпичная кладка; 5 — теплоизоляционная обмазка; 
р - диаметр наружной стенки кирпичной кладки; ©, №, №, №, У, У, У - геометрические 
размеры печи 


Температурное поле в печи должно соответствовать технологическим условиям обработки 


щебня. В связи с этим задача проектирования системы управления температурным полем печи 
приобретает особую важность. 
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Постановка задачи. Для проектирования системы управления температурным полем кольцевой 
роторной печи необходимо решить следующие задачи: 

- описать математическую модель температурных полей кольцевой роторной печи; 

- построить передаточную матрицу объекта управления, связывающую вектор входных 
воздействий (тепловые потоки форсунок) и вектор функций выхода (температуру в точках уста- 
новки датчиков); 

- рассмотреть возможные методы синтеза системы управления температурным полем рас- 
сматриваемого объекта; 

- синтезировать систему управления. 

Рассмотрим решение первых трех задач. 

При разработке дискретной модели температурных полей кольцевой роторной печи сде- 
лаем следующие допущения: 

- будем полагать, что температура на входе и выходе камеры 1 поддерживается посто- 
янной; датчики измерения температуры 3 расположены внутри печи в точках Х!-Х, У=У› , О„ 
(1=1.,... 5); 

- нижняя часть нагревательной камеры и боковые поверхности теплоизолированы; ще- 
бень рассыпан равномерно по дну печи; 

- скорость движения щебня в нагревательной камере 0,0222 м/мин; 

- управляющим (входным) воздействием служит тепловой поток, вырабатываемый путём 
сжигания природного газа в форсунках 2, а его плотность вычисляется по формуле О;=Р; /Хз У, 
Р— мощность, выделяемая /-й форсункой. 

Геометрические параметры рассматриваемой печи: 2=9,88 м; Х!=1,40 м; №=0,20 м; 
Х№=0,20 м; №=0,30 м; =0,4 м; %=0,15 м; =0,8 м; ©,=30°. 

Дискретная модель. Схема дискретизации математической модели [1] рассматриваемого объ- 
екта управления приведена на рис. 2. В рассматриваемой задаче будем полагать, что по выбран- 
ной координате шаг дискретизации остается постоянным. 
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Рис. 2. Схема дискретизации математической модели объекта управления: 
Т!—75- температурные поля различных сред; п — точки дискретизации по радиусу (0/2) (м=1...., из); 
Дг -— шаг дискретизации по радиусу г п -— точки дискретизации по координате у (1=1...., Пз); 
Ду - шаг дискретизации по у; у — точки дискретизации по © (у=1...., уд; ДО — шаг дискретизации по © 


Согласно приведённой схеме дискретизации, математическая модель объекта управления 
в дискретном виде записывается следующим соотношением [1]: 
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где Ат - шаг дискретизации по времени; а; — коэффициент температуропроводности Ёй среды, 
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а, =—^; с- теплоемкость; /.,- коэффициент теплопроводности; ‘у — плотность, 1<у<ук 
су 


1 


/- номер среды, /=1, 2,..., 5; интервалы изменения параметров |, п выбираются в соответ- 
ствии со схемой на рис.3 (например, для /=1; и1<Н<Н>, П.<П<п?). 
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Рис. 3. Схема дискретизации функций выхода 


Дискретные аналоги граничных условий [2] отражают неразрывность температур и тем- 


пературных полей на границах раздела сред. Для фазовой переменной Т; они могут быть записа- 
ны в виде следующих соотношений 
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ра<Н<р», 1<У<ук, 
где ^, — коэффициент теплопроводности первой среды. 


Аналогичными соотношениями записываются граничные условия и для других фазовых 
переменных: 


т == 0; а — 0; Н=1...., Из, П=1...., Пз. 


Дискретный аналог входного воздействия О, (/=1...., 5) на объект управления (см.рис.3) 
записывается в виде следующих соотношений: 


ее и Е Ра = О 
1 . 
Аг г 


П:<П<П}?, У; <У<Ул1; Е!..., 5, 
где О; -тепловой поток, выделяемый Рй форсункой (]=1,..., 5). 

Функцией выхода служат значения температурного поля в точках установки датчиков Ти 
па; у (напаруа, (1=1,..., 5) - координаты точек установки датчиков) (рис.3), О; — количество теп- 
лоты, выделяемое }-й форсункой (1=1,..., 5). 

При моделировании объекта были использованы следующие значения геометрических и 
физических параметров: 


р=9,7; Ж=1 42: №=0,22; №=0,4; 
у! =0,8; >=0,4; У3=0,12; У4=0,8; 
№1=3; №2=15; Из=17; 
П1=3; П2=14; Пз=18; 
А/=0,1; Лу=0,06; 9=2,5°;у,=120; 


у!=23, \=25, уз=41, у 4=43, \5=59, \=61, 
\у7=77, У в=79, у9=95, У10=97; 
4=14, па=8, (1..5); 
уа=8, уф=26, у@=44, уа=62, у%=90. 


0=10007 

а1=0,014 — воздух; 
а2=0,016 — кирпич; 
аз=0,008 — щебень; 
[41 =18 — воздух; 
142=23 — кирпич; 
143=20 — щебень; 


Г=0,00037 — —м/с. 
По результатам расчетов получена передаточная матрица И/ рассматриваемого объекта, 
связывающая вектор входных воздействий (тепловые потоки форсунок) и вектор функций выхода 
(температуру в точках установки датчиков): 






























































0,015 ск 0,0125 5 0,0095 к 0,0065 ив 0.004 пик | 
206051 502051 6900$ +1 3400$ +1 8600$ +1 
0,051 с 0,050 ж 0,0385 ив 0,0265 эк 0.0155 мы 
238051 44605 +1 69205 +1 3520$ +1 8940$ +1 
п’_| 0,042 зв 0,0765 45 0,068 5 0,0465 мв  0.027 виз | 
5140$+1 42765 +1 58805 +1 770051 8200$ +1 
0,0315 в _ 0,059 вк 0,083 5 0.067 5 0,039 зы 
74405 +1 6740$ +1 49845 +1 57005 +1 680051 
0,016 6-08 0,030 6798 0,042 ©7248 0,0555 6-99 0,056 675 
[8600$ +1 34805 +1 6920$ +1 4700$ +1 2486541 | 
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Выделим два сложившихся в настоящее время подхода для построения замкнутых систем 
управления рассматриваемым объектом, передаточная матрица которого приведена выше: реше- 
ние методами сосредоточенных систем и решение методами систем с распределёнными парамет- 
рами. 

Решение методами сосредоточенных систем основывается на свойстве диагональной до- 
минантности передаточных матриц объекта управления [2,3]. Положим, что в результате экспе- 
риментальных исследований получена передаточ- 
ная матрица, связывающая /Й вход с #м выходом: 


(5) =, ,(5)| Исследуем полученную переда- 


точную матрицу. Полагая 5=), где « — круговая 
частота, определим модули и фазы элементов мат- 
рицы И/. Если модули диагональных элементов мат- 
рицы И/ намного больше суммы модулей остальных 
элементов соответствующей строки, то матрица И/ 
обладает свойством диагональной доминантности [3]. 
Исследование диагональной доминантности может 
Рис.4. Построение спектра Гершгорина: быть осуществлено с помощью спектров Гершго- 
а; — заданные значения круговой частоты &;  рина для заданной передаточной матрицы. При 
АХ _ А, этом для каждого значения «= «„ могут быть по- 

т строены круги Гершгорина, объединяя которые, по- 

К= >. [(пк ИИ, (фо)? = Ве(И/, д)"; лучим спектр Гершгорина (спектр Гершгорина, по- 


ее в я строенный для Ай строки диагонально-доминантной 
| ]= .., т, тхт- мерность матрицы 
иен О ре матрицы, приведен на рис. 4). 





Для диагонально-доминантных матриц разработана частотная методика синтеза много- 
мерных систем управления [3]. В рассматриваемой задаче диагональная доминантность отражает 
взаимовлияние }-го входного воздействия на АЙ выход. 

На рис.5 показан спектр Гершгорина для первой строки передаточной матрицы рассмат- 
риваемого объекта (\/). Аналогично могут быть построены спектры Гершгорина и для остальных 
строк матрицы. 


т 





Рис. 5. Спектр Гершгорина для первой строки матрицы И/ 
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Выводы. В рассматриваемой статье описана дискретная модель температурных полей кольце- 
вой роторной печи и построена передаточная матрица рассматриваемого объекта. С помощью 
спектров Гершгорина исследованы характеристики передаточной матрицы. Результаты расчётов 
показали, что передаточная матрица рассматриваемого объекта не обладает свойством диаго- 
нальной доминантности, т.е. взаимосвязи между м входным воздействием и &й функцией выхода 
достаточно существенны, и их нельзя не учитывать. Следовательно, методика синтеза регулято- 
ров многомерных систем управления, опирающаяся на диагональную доминантность передаточ- 
ной матрицы объекта управления, не может быть использована. В рассматриваемом случае может 
быть использована методика синтеза систем с распределёнными параметрами, описанная в [4]. 
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